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ABSTRACT 

We present a spatially resolved comparison of the stellar-mass and total-mass surface distri- 
butions of n ine early- type gala xies. The galaxies are a subset of the Sloan Lens ACS survey 
(or SLACS: [Bolton et ani2006h . The total-mass distributions are obtained by exploring pixe- 
lated mass models that reproduce the lensed images. The stellar-mass distributions are derived 
from population-synthesis models fit to the photometry of the lensing galaxies. Uncertainties 
- mainly model degeneracies - are also computed. Stars can account for all the mass in the 
inner regions. A Salpeter IMF actually gives too much stellar mass in the inner regions and 
hence appears ruled out. Dark matter becomes significant by the half-light radius and becomes 
increasingly dominant at larger radii. The stellar and dark components are closely aligned, but 
the actual ellipticities are not correlated. Finally, we attempt to intuitively summarize the re- 
sults by rendering the density, stellar-vs-dark ratio, and uncertainties as false-colour maps. 

Key words: gravitational lensing — dark matter — galaxies: elliptical and lenticular, cD — 
galaxies: evolution — galaxies: haloes — galaxies: stellar content 



> 1 INTRODUCTION 

o\ 

In the current paradigm of galaxy formation, the building block of 
T-H structure is a dark matter halo, consisting mainly of non-baryonic 
dark matter together with ~ 15% baryons in the form of gas. 
^~>) , Dark halos are thought to originate from the collapse of primordial 
density fluctuations, growing unimpeded until virial equilibrium 
J — _ is reached. Within these halos the baryonic component dissipates 
. energy, collapsing further towards the center and eventually form- 
ic/ ' ing the visible galaxy. Subsequent mergers redistribute the matter 
. ^ ] within halos. 

• To work out the details of the basic picture, which goes back to 
$_( ' IWhite & R ees ( 1978), it is essential to determine the connection be- 
5^ , tween the visible galaxies and dark halos. Strong gravitational lens- 
ing by galaxies is potentially a very useful way of doing this, since 
the total mass of a lensing galaxy is relatively easy to constrain. 
Also, since lensing tends to be more effective for distant galax- 
ies (ziens ~ 0.1 to 1), it nicely complements the stellar-dynamical 
techniques applicable in nearby galaxies. 

The difficulty with lensing galaxies (that is to say, galaxies 
producing multiple images of background sources) is that they are 
relatively rare. Till recently, only about 80 were known. But re- 
cently, 28 new galaxy lens es have been disc overed by the Sloan 
Lens ACS Survey (SLACS: lBolton et alj2006l) thanks to a new sur- 
vey strategy, which eventually may double the number of strong 
lenses. The method is to select galaxies from the Sloan survey 

* Current address: Mullard Space Science Laboratory/UCL, Holmbury St. 
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jYork et alj|2000l) that have emission lines indicating high-redshift 
background objects, and then to image these candidates using the 
Advanced Camera for Surveys on the Hubble Space Telescope. 

A basic analysis of a sample of galaxy lenses is to fit sim- 
ple lens models and then comp ute M/L and its ev olution with 
redshift. This was first done by iKeeton et alj l ll998l) . An exten- 
sion is to place lensing galaxies on the Fundamental Plane of el- 
lipticals, using a measured or model-derived velocity dispersion 
(Kochanek et al. 2000; Rusin et al. 2003; Treu et al. 2006). Other 
work jTreu & KoopmansI 12004 ; Koopmans et al. 2006) compares 
lens models with the measured dispersions to constrain the mass 
profiles of the galaxies. These studies have found no unexpected 
features or trends with redshift, and argue in favor of passive evo- 
lution. 

A more detailed analysis involves modeling both the lens mass 
distribution and the stellar population. The star-formation history is 
not well-constrained by the observed fluxes and colours and must 
be marginalized over, but the stellar mass is fairly insensitive to 
model assumptions apart from the initial mass function (IMF). The 
lensing mass distribution, when aggressively modeled, turns out to 
have much larger uncertainties than simple models assume; never- 
theless, fliejinceitainties_can be estimated and useful conclusions 
drawn. iFerreras et alj ( 12005 ) found massive ellipticals to show a 
transition from no significant dark matter within Re to dark-matter 
dominance by 5-Re, whereas lower-mass galaxies showed no signif- 
icant dark halos even at 5Re. The radial gradient in the dark- matter 
fraction agreed w ith th e results qnjiearby galaxies derived from 
stellar dynamics jNapolitano et al.l200A . 

In this paper we extend the detailed comparison of stellar and 



© 2007 RAS 



2 /. Ferreras, P. Saha and S. Buries 



s 

o 



11.5 



11.0 



10.5 



10.0 



9.5 



- 
















1 ! ! ! ! 
























= 0.5 
























Z 
















































































































































































































































































































































































_ 




































_ 




































_ 




































_ 




































_ 








































z = 





1 




































































































































14W 


= 20. 


0. 






























i i 

























2 3 

(F435W-F814W), 




Figure 1. Stellar masses derived from the available photometry. The pre- 
dictions of the Bruzual & Chariot (2003) models are shown for a Chabrier 
(2003) IMF. A galaxy with apparent magnitude F814W = 20 is con- 
sidered at two different redshifts as labelled. The shaded regions span 
the model predictions when assuming a wide range of r models (i.e., an 
exponentially decaying star formation history at fixed metallicity). The 
dark/light shaded regions con'espond to a formation redshift of zp = 5 
and 2, respectively, and all span a range of star formation timescales 
(-1 < logT(Gyr) < +1) and metallicities (-1 < log Z/Zq < +0.3). 



total mass to two dimensions, using a subsample of the SLACS 
lenses. All the SLACS objects have a background galaxy that is 
lensed into two or four extended images. In nine of the objects, 
we can identify small features within the extended images. For 
lenses showing point-like multiply-imaged features there is a well- 
developed technique for reconstructing t he projecte d mass distri- 
bution, along with uncertainty estimates jSaha & Williams 2004 ). 
Accordingly we take these nine lenses as our sample. The stel- 
lar mass content is estimated by combining the available pho- 
tometry of the lensing ga laxies with population-synthesis models 
( iBruzual & Charlotll2003h . 



2 MAPPING THE TOTAL MASS 

We use the pixelated lens reconstruction method implemented in 
the PixeLe ns program^l The algori thm is motivated and described 
in detail in lSaha & Williamsl ( l2004h . but basically it consists of two 
ideas. 

The first idea is to express the lensing data as a set of lin- 
ear constraint equations on the mass distribution. We assume that 
the centroids of multiple images are measured with negligible er- 
ror, which — at HST resolution — is a very good approximation. 
Then the remaining unknowns are the source position and the lens 
mass distribution, and both of these enter th e lens equation linearly 
(cf. Equation 5 in lSaha & Williams! 1 1 997l) . The lens is expressed 



^ Available from jhttp : / /www . qgd ■ uzh . ch/pro jects/pixelens/| 



Figure 2. Projected mass maps of the lensing galaxies. The gray scale rep- 
resents Sgtcl while the contours show Etot ■ The latter are in multiphcative 
steps of 10" "' (like a magnitude scale). The star symbols mark multiply- 
imaged features of the background source. The bar in each panel indicates 
2 kpc. 



as a superposition of mass tiles or pixels, although a superposition 
of basis functions is also possible (Trotter et al. 2000). A prior on 
the lens can be expressed as linear inequalities. Specifically, we 
require the mass distribution to (i) be non-negative, (ii) be cen- 
trally concentrated, with the local density gradient pointing 45° 
away from the center, (iii) be inversion symmetric (i.e. symmet- 
ric with respect to rotation by 180°; optional), (iv) have no pixel 
more than twice the sum of its neighbors, except possibly the cen- 
tral pixel, and (v) the circularly averaged mass profile to be steeper 
than J?~°'^, where R is the projected radius. Item (v) is assumed 
(for galaxies) since stellar dyna mics or other methods never give 
p(r) as shallow as r~^'^ (e.g.. iBinnev et al]|l99lh lOerhard et al.l 
I2OOII) . These various equations and inequalities give an underde- 
termined linear system analogous to Schwarzschild's problem in 
stellar dynamics. 

The second idea (suggested by ^ Williams & Saha !'2000) is to 
sample the mass maps allowed by the data and prior, through a 
Monte-Carlo method. The result is an ensemble of lens models. 
Uncertainties on any parameter of the lens can be derived from the 
model ensemble in the usual way. The ensemble-average will auto- 
matically satisfy the data and prior constraints, since they are linear. 
This makes the ensemble-average model a good choice for repre- 
senting a single model. For a detailed illustration of the ensem ble 
method, see the Appendix in ISaha. WilUams & Ferreraj | |2007|) . 

For each lens we computed an ensemble of 200 mass maps of 
21x21 pixels each. To derive radial profiles we circularly-averaged 
the mass maps. This procedure naturally yields ensembles of radial 
profiles, hence an estimate of the uncertainty. 

In the lens models the enclosed mass is best constrained at 
projected radii similar to the images themselves (roughly speaking, 
around the Einstein radius), whereas at smaller and larger radii, 
the enclosed mass becomes progressively more uncertain. This is 
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Figure 3. Extended images in the lens J 1 205+491. Main panel. Extended 
images predicted by our model for this lens, for a circular source of size 
0.15" having a conical light profile. Inset. The observed image with the 
lensing galaxy substracted off. 



simply the well-known steepness degeneracy in lensing theory (see 
e.g., Saha 2000). As a result, the enclosed mass profile and its error- 
bars have a characteristic butterfly shape. 

To what extent velocity dispersions can break lensing de- 
generacies remains an open question. Studies using velocity dis- 
persions as constraints on lenses tend to be optimistic. However, 
stellar- dynamical degeneracies associated with kinem atics are also 
known jValluri et alj2004ICretton & Emsellernll2004) and it is not 
clear that these will be orthogonal to lensing degeneracies. Also, 
current methods for incorporating kinematics into lens models as- 
sume spherical symmetry. For these reasons we do not include 
kinematics in the lens models. 

We do not attempt to fit the extended images. Now, a model 
fitted to point-like images automatically provide a predictions for 
arcs that would be generate d by a conical l i ght p rofile (i.e. a 
circularly symmetric model; I Saha & Williamsl j200L) ) and these 
predicted arcs generally do yield a good approximation to ob- 
served arcs. Hence disregarding the extended images does not 
sacrifice as much information as might at first appear. Fitting 
a better source distribution is a straightforward linear image- 
reconstruction problem if the lens model is fixed. But using the 
source structure to improve the lens fit is much harder. Several 
papers propose schemes that iteratively fit lens and source mod- 
els. But th ey either assume a restricted parameterized form for the 
lens (e.g., Wallington. Naravan & KochaiieS 1994 Dye & Warred 



l2005l : iBrewer & Lewisll2004 iKoopmans et aljl2006h or allow the 
iKoopmans 1120051 : 



lens model to be free form iKoopmans II2005I : ISuvu & BlandfordI 
l2006h but need it to be close to a starting guess. 



3 MAPPING THE STELLAR MASS 

The galaxy images are Wide-Field Channel/ACS snapshots with 
an exposure time of 420s in each of the F435W [B) and F814W 
(J) filters. For each system we worked on the best-fit smooth lens- 
galaxy image that was originally subtracted from the actual image 
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Figure 4. A comparison of the ellipticity of the total and stellar mass distri- 
butions. The dots give the values of the measured ellipticities whereas the 
orientation of the arrows represent the misalignment between both distribu- 
tions; a vertical arrow means perfectly aligned. 



in order to discern faint lensing features. The lens-galaxy images 
were further corrected for Galactic reddening using the available 
dust maps (Schlegel et^l. 1998 ) and a standard model for extinc- 
tion by Galactic dust ( lFitzpatrickll999() . We then spatially rebinned 
the photometry of each galaxy into the same 21x21 pixelation used 
for the lens models of that galaxy. These pixels are ~ 0.1" to 0.2" 
across, depending on the system, amounting to ~ 0.1 to 1 kpc at 
the lensing galaxy. 

For the stellar-population analysis we apply the so-called r 
models, which assume an exponentially decaying star formation 
history at fixed metallicity. Three parameters describe each model: 
the formation epoch (i.e., when star-formation starts), the formation 
timescale and the metallicity. We explore a grid of 64 x 64 x 64 
models and allow in our estimates only those models whose pho- 
tometry is compatible with the observations, treating the 21 x 21 
pixels described above as independent. One colour is not sufficient 
for an estimate of the stellar ages and metallicities, but is enough 
to constrain the stellar mass content, because of the wide spectral 
range covered by B — I. Figure [T] illustrates this point. In this 
figure we consider a Ichabrieii ( l2003h Initial Mass Function with 
formation epoch zf = 2 or zf = 5, viewed from z = 0.1 or 
0.5. The stellar mass corresponds to a galaxy with apparent magni- 
tude F814W = 20 at the labelled redshifts. The shaded regions 
show the range of B — / as the metallicity varies in the range 

— 1 < Io^Z/Zq < +0.3 and the star formation timescale varies as 

— 1 < logr(Gyr) < +1. The uncertainty in ages and metallicities 
result in a ~ 0.2-0.3 dex uncertainty in the stellar mass content. 

Our sample comprises early-type galaxies. In these sys- 
tems the presence of an overall old and coeval stellar population 
jStanford et al.ll993) suggests star formation took place in a strong 
burst that consumed most of the available gas jFerreras &Silkl 
l2003h . Other baryonic components, hot and cold gas or dust, are 
known to contribute a small fraction to the net baryon budget in 
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Table 1. Summary of the galaxy sample. (Indicated errors are 90% confidence.) 



ID 


Rl Redshifti 


cr^ Rm 


Mtot(< Rm) 


Mstar(< Km) 


ba^Re/G 




kpc 


km/s / Re 


lOi°M0 


10"' A/0 


10"' M0 



J003753.21-094220.1 


6.99 


0, 


.195 


265 


1, 


.66 


52.3; 


J073728.45+321618.5 


10.10 


0, 


.322 


310 


0, 


.79 


38.3; 


J09 1205. 30+00290 1.1 


9.47 


0, 


.164 


313 


0, 


.95 


71.2; 


J095629.77+5 10006.6 


8.85 


0, 


.241 


299 


1, 


.20 


66.4; 


J120540.43+491029.3 


8.03 


0, 


.215 


235 


1, 


.04 


38.6; 


J133045.53-014841.6 


1.28 


0, 


.081 


178 


2, 


.19 


4.! 


J163602.6I+470729.5 


5.41 


0, 


.228 


221 


2, 


.01 


52.9; 


J230053. 14+002237.9 


6.43 


0, 


.229 


283 


1, 


.24 


40.9; 


J230321.72+142217.9 


8.13 


0, 


.155 


260 


1, 


.17 


49.8" 



+35.4 
-14.5 
+ 11.8 
-3.9 
+ 14.5 
-21.7 
+25.7 
-16.7 
+ 16.6 
-5.5 
n+3.5 
^-1.7 
+ 16.6 
-14.3 
+ 14.2 
-3.1 
+22.8 
-15.5 



44.3 
40.2 
44.7; 
41 
50.5 
3.8 
32.2 
28.0 
37.5 



+5.1 
4.3 
+4.1 
3.6 
+4.6 
4.0 
+4.9 
-4.0 
+5.1 
4.6 
+0.4 
0.4 
+3.8 
3.2 
+3.0 
2.6 
+4.0 
3.5 



57.1 
112.9 

107.9 
92.0 
51.6 
4.7 
30.7 
59.9 
63.9 



1 Data from Bolton et al. (2006). 



these systems jRobeits & Havneslll994h . Hence, we can take the 
stellar mass content in each modeling pixel as the baryon content. 



4 COMPARISON OF STELLAR AND TOTAL MASS 

With the methodology described in the previous two sections, we 
generated non-parametric 2D maps of the total and stellar surface 
density on the same pixelation. Figure |2] shows both distributions, 
the gray scale being Estoi derived from population synthesis and 
the contour maps indicating the Etot derived from lensing. (Both 
maps have uncertainties, shown in later figures, but not here.) Seven 
of the lenses have the inversion-symmetry constraint mentioned in 
Section[2] For the two lenses where the images are well distributed 
in position-angle, J0912 and J1636, we did not impose inversion 
symmetry and the models are allowed to be lopsided. The asym- 
metry in these two lenses is at the level of 10-15% which given the 
uncertainties shown later in Figure[5]is probably not significant. 

We noted in Section [2] that a model fitted to point-like fea- 
tures can be used to predict extended images to some degree. 
For one case this can be achieved with a simple piece of com- 
puter graphics. If one draws the arrival-time contours of the point 
source with a very close contour spacing, the resulting pattern mod- 
els the lensed image of an ext ended source with a conical light 
profile dSaha & Williamsll200lh . Figure [3] illustrates for one lens, 
J1205+491, which occupies the central panel of Figure|2] The white 
curves are the contour lines, and two minima and two saddle points 
are discernible; these are the locations of the point-like images. In 
most of the figure the contour lines are so close that the result- 
ing pattern is almost completely white, and has been airbrushed 
out. In the region shown, however, the white-on-black pattern 
closely resembles the observed arcs. The white lines have width 
w = 0.01 arcsec on the scale of the figure. The contour spacing in 
scaled units (see e.g.. Equation 2.3 in iBlandford & Naravaij[T983) 
has been chosen as Ar = 0.0015 arcsec^. Th e implied source 
radius in this type of model is Ar/w or 0.15". iKoopmans et al.l 
( l2006h have much more detailed source maps (but much simpler 
lens models) and this source size is typical. The main conclusion 
from Figure [3] is that most of the lensing information in the ex- 
tended images is already present in the point-like features. 

Returning now to examine Figure[2]again, comparing the stel- 
lar and total mass profiles suggests that 

(i) Estei and Etot tend to have aligned ellipticity, and 

(ii) Estoi falls off more steeply than Etot • 



To test the statement (1) above, we compute the average el- 
lipticity and position angle of Estoi and Etot using the second or- 
der moments of the surface density. Figure |4] compares these. The 
orientations appear to be almost perfectly aligned. The tendency 
of lens models to be oriented with the light is well known (e.g., 
lKeetonetal1ll997h . At the same time, the magnitudes of the ellip- 
ticities is not correlated; in some galaxies the stellar part is rounder, 
in some galaxies the dark matter is rounder. 

Statement (ii) above is addressed by Figure [S] which shows 
the circularly-averaged enclosed mass profiles, with 90% Bayesian 
confidence intervals. The error bars in the total mass derive from 
the ensemble of lens models, and have a characteristic butter- 
fly (or bow-tie) shaped envelope as explained in Section 2. One 
envelope of the butterfly shape is expected to be approximately 
Mtot{R) <x i?^"'', because the prior requires Etot to be steeper than 
T?""'* (see §2). This is indeed the case, as illustrated in the bottom- 
middle panel of the figure. The other envelope has a more com- 
plicated origin, having to do with the steepest non-negative profile 
that can reproduce the observed image positions. The uncertainty 
in the stellar mass arises from the variety of metallicities and star- 
formation h istories compat ible with the available photometry (see 
Figure[T). A lChabriej 12003) form for the IMF is assumed, with the 
effect of changing to a Salpeter IMF shown for one galaxy — we 
will return to this question below. 

The strongest inference from Figure [5] is that dark matter is 
located in halos. While there is diversity in the details, the gen- 
eral pattern is that stellar mass dominates in the inner regions with 
the dark-matter fraction consistently increasing with radius0This 
conclusion is as expected based on our knowledge of stellar and 
gas dynamics in nearby galaxies, but is emerging here from a com- 
pletely different method and data set. 

Our sample spans a similar range of masses and velocity dis- 
persion (see Table[Tll. Except for J1330, the total mass extends only 
over a factor of two. All galaxies appear baryon dominated at their 
centers, whereas outside the half-light radius the profiles show a 
wide range of distributions. Most of the galaxies show a signifi- 
cant contribution from dark matter. Interestingly, the only galaxy 



^ We remark that the above result cannot be an artefact of the lens- 
modelling prior because (a) the prior has no influence on the normalization 
of the lensing mass, which depends only on image positions, redshifts, and 
cosmological parameters, and (b) a spurious halo would amount to a Stot 
profile that was not steep enough, whereas the prior imposes a minimum 
steepness. 
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Figure 5. Cumulative mass profiles of the lensing galaxies. Black circles 
and grey triangles denote total and stellar mass, respectively. The error bars 
show the 90% confidence region. The effective radius {Re ; measured from 
the optical image) is given as an arrow for each galaxy. The enclosed mass 
is best constrained in the region of the images, with the error bars enlarging 
at smaller and larger radii. Outside of i?e all the galaxies have a significant 
contribution from dark matter although some systems, most notably J2300, 
present a very high dark matter contribution. A Chabrier IMF is used for the 
stellar mass estimates. The grey crosses on the upper left panel correspond 
to the stellar mass for a Salpeter IMF, which gives an unphysical Stot < 
Sstcl in the inner parts of the galaxy. The dashed curve in the bottom- 
middle panel shows R^-^ , the steepest enclosed-mass profile (equivalently 
shallowest density profile) allowed in our models. The normalization of the 
dashed curve is arbitrary. 



with no evidence of dark matter to the radius observed (J0737) is 
also the one where the observations are all interior of the half-light 
radius. This suggests that dark halos become significant roughly 
around the half-light radius. 

As mentioned above, we have assumed a Chabrier IMF for the 
stellar population, and other realistic choices of IMF differ from 
this one by a small fac tor not includ ed in the error bars. The simple 
power law defined bv lSalpeteil ('l955) and traditionally used as an 
approximation to the IMF gives stellar masses ~ 50% higher. For 
some galaxies presented here, the implied stellar mass exceeds the 
total mass as extracted from lensing. The upper-left panel of fig- 
ure |5] shows the stellar mass profile of J0037 for a Salpeter IMF, 
which is clearly incompatible with the lensing results. Hence, we 
can infer from our results that a Salpeter IMF is too bot tom-heavy. 

Our conclusion agrees with lCappellari et^ ( l2006h . who used 
stellar dynamics to estimate the total mass in a sample of 25 E/SO 
galaxies. 

We can compare the aperture masses (stellar and total) with 
previous work. The aperture radius Rm we take as 2i?max — R-mi-a, 
where 7?max and i?min are the projected radii of the outermost and 
innermost lenses images. This Rm approximates the radi us within 
which lensing usefully constraints on the mass. In Fer reras et al.l 
( l2005l) the aperture radius was set for each galaxy by seeing where 
the error bars became too large to be useful; such a procedure 
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Figure 6. Stellar and total mass compared in lensing galaxies and in nearby 
galaxies with stellar dynamics. Main panel. Dots show the aperture masses 
A/stcl I Aftot . (The aperture is the region where the total mass is usefully 
constrained by lensing.) Black dots come from the present work using 
SLACS lenses, and gray dots from a similar analysis on a sample of CAS- 
TLES galaxies jperreras et alJllOO^ . The stars are measurements based 
on dynamical models applied to (non-lensing) elliptical galaxies observed 
with the SAURON Integral Field Unit (Cappellari et al. 2006). The dashed 
line corresponds to A/tot = A'/gtoi, while the dashed line corresponds to 
A/tot ^^s\ol' would give the observed tilt in the Fundamental 

Plane. Top panel. The aperture radius Rm relative to the half-light radius 
R&. 



tends to give slightly smaller aperture radii, but we will disregard 
the difference here. Figure |6] shows the aperture Mtot and Mstd 
in this work alon g with the earlier lensing work (labelled CAS- 
TLES: iFerreras e t al. 2005|}_ and the stellar-dynamical work (la- 
belled SAURON: ICappellari et al.ll2006 ). The galaxies from these 
three independent data sets all show a clear trend towards a larger 
amount of dark matter in more massive galaxies. The observed tilt 
of the Fund amental Plane suggests a scaling of Mtot oc Afg\,?i 
jFerreras & Silk 2000, 2003), which appears consistent with all the 
data. However, the present sample is too small to test a scaling 
law, especially because the range of aperture masses is small. Also 
noticeable in the figure is that the SLACS galaxies tend to have 
Rm / Re — 1, unlike the CASTLES galaxies which cover a larger 
range. This is a selection bias due to the SLACS survey strategy. 

In Table [T] we include an estimate of the dyn amical mass, us- 
ing th e popular relation Mvir = ^ReO'^/G (e.g.. lCappellari et al.l 
I2OO6I) . In the table Rm represents, as explained above, the aperture 
radius usefully mapped by lensing. Notice that for J0737 {Rm = 
0.8i?e) the dynamical mass estimate is much larger than the lens- 
ing or stellar masses within Rm, whereas J1330 (Rm = 2.2i?e) 
is better mapped by the lensed images and give compatible results 
between Mtot and the dynamical mass. 

Finally, Figure [7] shows stellar and dark-matter maps in false 
colour. Here the brightness represents Egtci; the colour represents 
the dark-matter fraction (red means all stellar, blue means all dark); 
the hue represents the uncertainty (pale means more uncertain). The 
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Figure 7. A false colour map of stellar and total mass in galaxies. The surface mass density is show by intensity (black: low 
density; bright: high density). Dark matter fraction is represented by colour (blue: dark matter dominated; red: stellar matter 
dominated). Uncertainty is encoded by hue (pale: high uncertainty, bright: low uncertainty). The colour bar on the left shows 
a range of mass ratios for two different uncertainties in the measurement of Sgtel/^tot: 0% (exact measurement; left) and 
20% en'or (right). The horizontal bar in each panel indicates 2 kpc. 



Appendix explains the colour-coding scheme more precisely. At 
the centers of the galaxies we see bright white, meaning high den- 
sity but with large uncertainty. This turns into pale red, indicating 
mainly stellar mass. Further out we see faint pale blue, indicating 
mainly dark matter. This figure, though admittedly ad hoc in its 
choice of colour coding, provides for the first time an intuitive vi- 
sual image of the distribution of dark matter in individual galaxies. 



5 CONCLUSIONS 

The potential usefulness of early-type lensing galaxies for under- 
standing the interdependence of baryons and dark matter in galaxy 
formation and evolution is widely appreciated. Lensing can be used 
to map the total mass and starlight can be used to map the baryons. 



In this pa per we do this for nine galaxies from the SLAGS sur- 
vey by Bolt on et al. H2OO6). From the lensing data we derived free- 
form pixelated models for the total mass, and from the galaxy pho- 
tometry we computed stellar population-synthesis models. In both 
casses we generated Monte-Carlo ensembles of models, in order to 
marginalize over unknowns such as lensing degeneracies and star- 
formation histories, thus obtaini ng realistic uncertaint ies. The tech- 
nique is basically the same as in lFerreras et alj jZOOSi) . but whereas 
the earlier work only compared radial profiles now we compare 
stellar and total mass in 2D. Related work has been done on larger 
samples of galaxies, but is limited to fitting simple parameterized 
models for the lens, and does not model the stellar populations at 
all. 

The 2D mass maps are shown in Figures |2] and |7] The former 
overlays contours of Etot on a grayscale of Sstci • The latter shows 
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the same information with uncertainties as well, all encoded in false 
colour: red for stellar mass, blue for dark, and pale versus coloured 
for uncertainty. It is evident that (a) these galaxies are dominated 
by stars in the inner regions, but mainly dark matter in the outer 
regions, and (b) stellar and dark components are well-aligned, but 
neither has a simple elliptical shape. 

One can of course still compute an ellipticity defined as a mo- 
ment, and this is shown in Figure |4] We see that ellipticities of the 
stellar and total mass are uncorrelated in magnitude but are almost 
perfectly aligned. It would be interesting to see if this is true of 
galaxy-formation simulations. 

The profiles (Figure O show that dark matter halos begin 
to dominate around the half-light radius, although some galaxies 
seem to present a stronger contribution from dark matter even in- 
side i?e (e.g., J2300). The present sample is too small to extract 
any strong correlation of the dark matter distribution with global 
properties such a s total mass or lumi nosity. However, the trend 
found previously jFerreras et alj[20o3) , namely that there should 
be more dark matter in more massive galaxies, with a scaling of 
roughly A/tot oc Ml^^i (which is equivalent to the tilt of the Fun- 
damental Plane) is compatible with the combined data of lensing 
galaxies (labelled CASTLES and SLACS) along with the dynami- 
cal analysis of local ga laxies with the SAURON integral field unit 
jCappellari e"tai]|2006h . 

The top-left panel of Figure [5] also shows that a Salpeter 
(1955) IMF cannot be used to estimate stellar masses as the pop- 
ulation synthesis models predict too much stellar mass compared 
to the total mass obtained from our lensing studies. This result also 
agrees with the analysis of the SAURON sample. 

We emphasize that stellar and total masses are obtained from 
different data through models of different physical processes. There 
is no tuning to give similar results. That the stellar and total densi- 
ties come out compatible at the centers — where the baryon content 
is expected to dominate the mass budget — indicates that systematic 
effects are not significant. 
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APPENDIX A: FALSE-COLOUR MAPS 

For each galaxy we have sky-projected maps of three quantities: to- 
tal mass, uncertainty, and stellar-mass fraction. How can we some- 
how encode these into red, green, and blue, thus making a false- 
colour map? 

In general, suppose we want to represent the following three 
quantities: 

A: an amplitude 

/: a ratio of components (Al) 
A: a fractional uncertainty 

where / and A both vary between and I. A plausible mapping 
into intensities r,g,b of red, green, and blue is 

r/A = iA+(l-A)/ 

g/A = iA (A2) 
b/A = iA+(l-A)(l-/). 

This gives total intensity r + g + b — A regardless of / and A. 
If A = 0, the colour will vary from blue at / = to magenta at 
/ = 0.5 to red at / = 1. As the uncertainty A increases, these 
colours will become paler, turning into white at A = 1. 

The above scheme on its own is, however, not enough to pro- 
duce useful false-colour maps. That is because the response of eyes 
to colour is highly nonlinear, contrary to what Equation l lA2b pre- 
supposes. In practical image-processing, heuristic scale stretchings 
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are always necessary. After some experimentation, we found the 
stretchings 

A = ^li?, f=[^f\ A=(^)'''.(A3) 

to be useful. Estci and Etot represent the surface mass density from 
the photometric and lensing analysis, respectively, in a given pixel. 
A is the fractional uncertainty in the total surface mass density, 
which dominates the error budget. Figure[7]shows the results. 



